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Electric dipole moments of diamagnetic atoms of experimental interest are calculated using the 
relativistic Hartree-Fock and random-phase approximation methods, the many-body perturbation 
theory and configuration interaction technique. We consider P,T-odd interactions which give rise 
to atomic electric dipole moment in the second order of the perturbation theory. These include 
nuclear Schiff moment, P,T-odd electron-nucleon interaction and electron electric dipole moment. 
Interpretation of a new experimental constraint of a permanent electric dipole moment of ^^^Hg 
[W. C. Griffith et al, Phys. Rev. Lett. 102, 101601 (2009)] is discussed. 
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I. INTRODUCTION 

Search for a permanent electric dipole moment (EDM) 
of particles, violating both parity and time-reversal in- 
variance, has a long history (see, e.g., [l|). The standard 
model predicts tiny EDMs which cannot be detected at 
the present level of experimental accuracy. However, dif- 
ferent extensions of the standard model (such as, e.g., 
supersymmetry) predict much larger EDMs of the parti- 
cles that, in principle, could be found using the modern 
experimental technique. A reveal of such EDMs would 
unambiguously lead to a contradiction with the standard 
model. 

According to 0, 0] previous limits on EDM produced 
stringent constraints on electroweak bariogenesis and 
models of CP-violation, but fall short of ruling out the 
simplest generic extensions of the standard model. It 
was stated that the next generation of EDM experiments 
should be sufficiently sensitive to provide a conclusive 
test. 

A very significant step ahead in this direction has 
been done in a recent work The authors reported 
new result obtained for a permanent EDM of ^^^Hg to 
be d(i99Hg) = (0.49± 1.29stat ± 0.76,y,t) x IQ-^^ e cm. 
Though the EDM is still consistent with zero, its limit 
|(i(^^^Hg)| < 3.1 X 10~^^ |e| cm is an improvement of the 
previous Hg limit by a factor of 7. Motivated by this re- 
sult and implying future experimental progress we have 
performed calculations of different contributions to the 
atomic EDMs of ^^^Xe, ^^^Yb, ^^^Hg, ^iiRn, and ^^^Ra. 
For these atoms the experiments searching for EDMs are 
underway. 

The paper is organized as follows. In Sec.|TT]we discuss 
different types of P,T-odd interactions that can lead to 
an appearance of a permanent atomic EDM. In Sec. IIIII 
we describe the methods of calculations of the EDMs. 
We start our calculations from the relativistic Hartree- 
Fock method. Than we include many-body corrections 
using two different methods. First, we apply a simple 



random-phase approximation (RPA) for the closed-shell 
atoms. Second, we apply the configuration interaction 
(CI) combined with the many-body perturbation theory 
(MBPT) approach to valence electrons while use the 
RPA approach for the core. Sec. IIVI is devoted to an 
analysis and discussion of the results. We present the 
results obtained for different contributions to the atomic 
EDMs and compare them with other available data. In 
Sec. |V]we discuss the neutron and proton contributions 
to the total nuclear spin using the spherical shell model 
of a nucleus. Sec. IVII contains concluding remarks and 
two final Tables where the recommended values of the 
contributions to the EDMs of ^^^Xe, ^^^Yb, ^^^Hg, ^^Rn, 
and ^■^^Ra are gathered and the limits on CP-violating 
parameters based on new experimental limit for ^^^g 
are presented. 



II. GENERAL FORMALISM 

Our goal is to find the atomic EDM dat defined as 
dat = dat(F/F), where F = J -f I with J being total 
angular momentum. In this work we deal with the atoms 
with closed shells in their ground states. In this case 
J = 0, F = I, and dat = dat (I//). 

We consider several types of P,T-odd interactions be- 
tween particles leading to an appearance of an atomic 
EDM. We restrict ourselves to contributions to the EDM 
which occur in the second order of the perturbation the- 
ory. The EDM induced in an atomic state |0) due to an 
admixture of opposite-parity states that appears in the 
second order of the perturbation theory can be written 
as 



dat = 2 ^ 



K 



mK){K\H\0) 
Eo — Ek 



(1) 



where D = — |e| r is the electric dipole operator, e is the 
electron charge, and Ei are the energies of the states. 
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We will consider below: 1) the tensor-pseudotensor 
P,T-odd electron-nucleon (e-N) interaction, 2) the 
pseudoscalar-scalar P,T-odd e-N interaction, 3) the nu- 
clear SchifF moment, and 4) the interaction of the elec- 
tron EDM with the internal nuclear magnetic field of the 
atom. 

It is worth noting that the operators describing all 
these interactions have certain similar features. All of 
them 1) are proportional to the nuclear spin I and 2) have 
strong singularity. The similar nature of these operators 
leads to a suggestion that the calculations depend very 
little on their specific form. It allows us to expect that 
the relative contributions of different many-body parts of 
calculations will remain approximately the same for all 
operators. As we will show below this suggestion is fully 
justified. 



A. Electron-nucleon P,T-odd interactions 

We start our consideration with brief reminder of the 
main features of the e-iV P,T-odd interaction leading to 
appearance of atomic electric dipole moments in the sec- 
ond order of the perturbation theory. A detailed descrip- 
tion can be found elsewhere 1 , [E] ■ This interaction has 
the following form (see, e.g. 



H 



G_ 



N 



+ Nij^Neel 



(2) 



Here and Cp are dimensionless coupling constant 
characterizing tensor-pseudotensor and pseudoscalar- 
scalar P,T-odd electron-nucleon interactions for the nu- 
cleon N; a^^ = {li^ilv - 7i^7p)/2, 75 and 7 are the Dirac 

matrices: 75 = ( 'j' \] and 7 cr 



^ -1 y ' \-cT Q 
In accordance with Eq. ^ we can represent the Hamil- 
tonian H as H = Ht + Hp, where in the coordinate 
representation (atomic units Ti — rrig = |e| = 1 are used 
throughout) 



Ht = iV2GCti{(tn) p{r) 



(3) 



G 1 

Hp = -= Cp 70 Vp{r) (cTjv). (4) 

y 2 zrUpC 



Here G is the Fermi constant, c is the speed of light (in 
atomic units c — l/a « 137) and rrip is the nuclcon mass. 
We denote 

Ct{o'n) = ( Ct E cTp + Ct E Cji ) , 



Cp{aN) = (c^^(Tp + C?^. 



where (...) means averaging over the nuclear state with 
the nuclear spin I. 

In Eq. (HI) we keep only the term in the lowest nonvan- 
ishing approximation in m~^. Note that in this equation 
the operator V acts only to p{r), where p{r) is the nu- 
clear density distribution. 

Since we are dealing with very singular operators the 
model of the nuclear density distribution can be impor- 
tant. To check this point we have carried out calculations 
for two models. In one of them the nucleus was treated 
as a charged sphere with the radius R, i.e.. 



4ttR 



B{R-r). 



In other model it was used the Fermi distribution 



p{r) = 



l + exp^:^- 



(5) 



(6) 



where po is the normalization parameter determined by 
J pdV = Z. We have found that the results obtained for 
each of these models were numerically very close to each 
other. 

Similar to the expression for H it is convenient to rep- 
resent dat as a sum of two terms dat = dat -I- d^ , where 
dat and d^ correspond to the operators Ht and Hp 
given by Eqs. Q and (jl]). Explicit expressions for dat 
and d^t can be derived from Eq. ([T]) by replacing the 
operator H to Ht and Hp, correspondingly. 

It is convenient to determine the quantities dat and 



as follows: dat^ = dt^^^ {(Tn) ^ d^^^ I/I. The coefficient 
of proportionality in this expression depends on a model 
of the nucleus. An accurate treatment of the nuclear 
structure is beyond the topic of this work. For a spherical 
shell model of the nucleus this coefhcient can be easily 
found for different atoms. We will discuss this problem 
in more detail in Sec. fVl 



jT,P 



jT.P 



B. The nuclear SchifF moment 

The Schiff moment is a nuclear moment violating both 
parity and time-reversal invariance. It is caused by P, T- 
odd nuclear forces and it takes into account screening of 
external electric field by atomic electrons. We use the 
form of the Hamiltonian for the interaction of atomic 
electrons with the nuclear Schiff moment suggested in 
Ref. i: 



3Sr , , 

HsM = — '^py'^i^ 



(7) 



where B = j p{r)r^dr and S is the Schiff moment vec- 
tor defined as S = S'(I//) with S being the coupling 
constant. 

Its contribution to the EDM of the closed-shell atoms 
has been considered in detail in 0, Q . As it follows from 
Q the results obtained in the frame of multiparticle ap- 
proach combining the CI with the MBPT agreed (within 
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10%) with the results obtained by the RPA method. Here 
we perform similar calculations mostly for consistency 
test and for completeness. 



C. The electron EDM 

An interaction of the electron EDM de with the elec- 
tromagnetic field strength f^i^ can be written in the rel- 
ativistically covariant form as 



(8) 



Here ip = V'Wo and is determined in Appendix 1X1 

Again we will consider here only effect appearing in the 
second order of the perturbation theory. It characterized 
by an interaction of the electron EDM with the magnetic 
field B created by the nuclear magnetic moment. The 
operator of this interaction {Hb) can be written as 



The magnetic field B can be represented by 
(M X r) 



B = V X 



(9) 



(10) 



3(Mn)n- M 



{r-R) + ^e{R-r), 



where 9{x) 



1, X > 
0, X < 

to the atomic EDM, can be written as d^j = 
found from Eq. ^ by replacing H Hb- 



and n = r/r. The contribution 



dg{I/I) and 



III. METHODS OF CALCULATION 

A. RPA for the closed shells 

Here we describe a simple method suitable for calcula- 
tions of atomic EDM for atoms with closed shell. On the 
first stage we solve Dirac-Hartree-Fock (DHF) equations 
in the approximation (i.e., including all electrons 
forming the ground state of the atom in a self-consistency 
procedure). 



(11) 



Here Hq is the relativistic Hartree-Fock Hamiltonian and 
■00 and Ec are single-electron wave functions and energies. 

At the next step we construct virtual orbitals. Dif- 
ferent techniques can be used for this procedure. One 
approach is to multiply the previous orbital of the same 
partial wave to a smooth function of r with subsequent 
orthogonalization of this orbital to all the rest orbitals. 
This method was described in detail in Refs. [ol. [l0|. 

Another method is to construct a basis set using the 
B-spline technique developed at the University of Notre 



Dame [llll. We use 60 B-splines of order 9 in a cavity 
of radius i?max = 30 as, where at is the Bohr radius. 
This relatively large number of B-splines is needed due to 
singularity of the P,T-odd operators. This requires very 
detailed description of the wave functions in the vicinity 
of the nucleus. 

Further, we consider an atom in external field and solve 
the RPA equations (self-consistent DHF in an external 
field) 



(12) 



where F is the operator of external field, 6VP is the cor- 
rection to the self-consistent potential due to the effect 
of external field. Index c numerates single-electron func- 
tions (ipc) of the closed-shell core. The RPA equations 
Eq. (fT2l) are solved self-consistently for all states in the 
core for all external fields involved in the problem. 



B. CI+MBPT 

A more sophisticated and accurate way to calculate 
atomic EDMs is to use the configuration interaction tech- 
nique for valence electrons while still use the RPA ap- 
proach for the core. This would allow us to check the 
accuracy of the RPA calculations. It is especially impor- 
tant for the atoms having two external s electrons: Yb, 
Hg, and Ra. The effect of external electrons on different 
properties of these atoms is large and accurate treatment 
of interaction between them is needed. 

We consider Yb, Hg, and Ra as atoms with two 
valence electrons above closed-shell cores [Is, 4/^''], 
[Is, 5d^°], and [Is, 6p^], respectively. In this paper 
we follow approach suggested in which combines the 
many-body perturbation theory with the configuration 
interaction method. We will refer to it as the CI-I-MBPT 
formalism. The MBPT is used to include excitations 
from the core into the effective Hamiltonian for valence 
electrons. After that the multiparticle relativistic equa- 
tion for valence electrons is solved within the CI frame- 
work to find the wave functions and the low-lying energy 
levels. 

In the Cl-f MBPT method, the energies and wave func- 
tions are determined from the time- independent equation 

H,s{En)^r. = EA, (13) 
where the effective Hamiltonian is defined as 



(14) 



Here iJpc is the Hamiltonian in the frozen core approxi- 
mation and S is the energy-dependent correction, which 
takes into account virtual core excitations. The operator 
S completely accounts for the second-order perturbation 
theory over residual Coulomb interaction. 

Since we are interested in calculating the atomic 
EDMs, we need to construct the corresponding effective 
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operators for valence electrons [l3|, ll5| . To do that we 
can extend the concept of the effective Hamiltonian H^g 
to other operators such as the effective dressed electric- 
dipole operator Doff and the P,T-odd operators. These 
operators account for the core- valence correlations. As 
in pure RPA approach of Sec. IIII A) we solve the RPA 
equations summing a certain sequence of many-body di- 
agrams to all orders of MBPT 13, 16, 17]. Since require- 
ments to the accuracy of calculations are not very high 
we disregard in this consideration small corrections like 
normalization and structural radiation. 

We perform the calculations in the CI-I-MBPT method 
in V^~^ and V'^~^ potentials. The former is a bit more 
"natural" for the lowest-lying odd-parity states of the 
considered atoms such as nsnp^'^P" (n = 6 for Yb and 
Hg, and n = 7 for Ra), because 6^1/2,3/2 (or 7pi/2,3/2) 
orbitals are constructed at the stage of solving DHF equa- 
tions for the configuration [core] nsnp. The latter is some- 
what simpler (e.g, due to an absence of the subtraction 
diagrams). We have checked that the final results in both 
potetials are in good agreement with each other. For this 
reason when we discuss results, we do not distinguish be- 
tween these potentials. 

Two different basis sets described in section UlI Al were 
used in the V^~^ and V^~^ approximations. First set 
[Tol ] was used in the V^^^ and B-spline set [ll| was used 
in the ^^-^ calculations. At the CI stage of the V'^'^ 
calculations the one-electron basis set for Yb included Is- 
13s, 2p-13p, 3d-12d, 4/-11/, and 5g-7g orbitals, where 
the core- and 5d and 6p orbitals are Dirac-Hartree-Fock 
ones and all the rest are the virtual orbitals. For Hg and 
Ra the basis sets were insignificantly larger. In all cases 
the basis sets were numerically complete and the full CI 
was made for two valence electrons. At the stage of the 
MBPT calculations a more extended basis sets, including 
more basis functions, were used. For instance, for Yb it 
included ls-26s, 2p-26p, 3d-25d, Af-Uf, and 5g-12g 
orbitals. 

All B-splines up to Imax = 5 were used to calculate S 
in the V^'"^ approximation. 18 lowest basis states above 
the core in each parial wave up to Imax = 3 were used on 
the CI stage of the V^~'^ calculations. 



IV. RESULTS AND DISCUSSION 
A. RPA for the closed shells 

In the frame of the RPA method discussed in Sec. IIIII 
we have performed the calculations of different contribu- 
tions to the EDMs of the diamagnetic atoms presented 
in Table [H For these atoms the experiments searching 
for the EDMs are planned or are underway. 

We can rewrite Eq. ([1]) as follows 



dat = 2^ 



a.k 



(fc|r|a)(fc|g|a) 



TABLE I: The isotopes of the atoms considered in this work. 
fj, are the magnetic moments expressed in nuclear magnetons. 





Z 


A 


I 




Xe 


54 


129 


1/2 


-0.7778 


Yb 


70 


171 


1/2 


0.4119 


Hg 


80 


199 


1/2 


0.5059 


Rn 


86 


211 


1/2 


0.60 


Ra 


88 


225 


1/2 


-0.734 



where the summation is over the quantum numbers of 
the one-electron core states "a" and excited states "k", 
and Ei are the one-electron energies. 



1. P,T-odd e-N interaction 

We start the discussion from tensor-pseudotensor and 
pseudoscalar-scalar e-iV P,T-odd interactions. Using the 
Wigner-Eckart theorem and going over to the reduced 
matrix elements (ME) we obtain for the contributions 
and dfj to atomic EDM 



2V2G^ ^ (A:|H|a)(fc||z7p(r)||a) 



a,k 



£k 



GCp 
3\/2 m„ 



a,k 



(fc|| r||a)(fc||n(dp/dr)7o||a) 



(16) 



(17) 



(15) 



The explicit expressions for the reduced MEs of the P,T- 
odd operators in Eqs. ([16]) and (fT7|) are given in Ap- 
pendix 13 

In Table HIl we present the values of c?^^. obtained for all 
isotopes listed in Table[T]in pure DHF approximation and 
including RPA corrections. Note that the RPA correc- 
tions must be included only for one operator in Eq. (jl6p 
or Eq. (|17p [l^. In other words, when we include the 
RPA corrections for the electric dipole operator r we 
must not include them for the P,T-odd operator and vice 
versa. Certainly both approaches should lead to the same 
result. It allows us to test consistency of the calculations. 

Our results for Xe are in good agreement with other 
data. Our RPA value for ^^^Hg is in excellent agreement 
with similar calculations by Martensson-Pendrill [l^ 
while the result obtained in Ref. [19| is somewhat larger. 
There is also a reasonable agreement between the result 
found in this work and the estimate obtained by Sushkov, 
Flambaum and Khriplovich from the analytically derived 
formula [13]. To the best of our knowledge there are no 
other data for Yb, Ra, and Rn to compare with. 

As is seen from Table [Hi inclusion of the RPA cor- 
rections leads to increasing atomic EDM. For the noble 
gases (Xe and Rn) the RPA corrections contribute at the 
level of 30%, while for the atomic Hg, Yb, and Ra, which 
have two s electrons above closed shells, the RPA correc- 
tions are much larger. In fact, they increase the EDM 
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TABLE II: The values of dft in units (10"2°Ct {o-]v)|e| cm) 
obtained in DHF and RPA approximations are presented. 
The results are compared with other data. 





i29Xe 


171 Yb 


199Hg 


2"Rn 225j^^ 


This work (DHF) 


0.45 


-0.70 


-2.4 


4.6 -3.5 


Ref.|jJ (DHF) 


0.41 




-2.0 




Ref.[21J (DHF) 


0.41 








Ref.[20| 


0.6 




-3.9 




This work (RPA) 


0.57 


-3.4 


-5.9 


5.6 -17 


Ref. [18] (RPA) 


0.52 




-6.0 




Ref.flQ] (RPA) 






-6.75 





2.5 times for Hg and 5 times for Yb and Ra as compared 
to the DHF results. The reason for this increase is that 
the two s electrons are loosely bound and can be easily 
excited. As a result, an account for the higher orders of 
the perturbation theory (like the RPA corrections) leads 
to significant change of the "bare" results obtained in the 
DHF approximation. 

The results obtained for df^ are listed in Table IIIII In 
Ref. [l^l Flambaum and Khriplovich suggested a method 
to establish the correspondence between contributions of 
the tensor-pseudotensor and pseudoscalar-scalar P,T-odd 
operators using the expressions for the reduced MEs of 
these operators (see Appendix [X)) . 

This correspondence can be obtained using the prop- 
erties of the wave functions /s, (?s, fpl/2^ ^^'^ 9pi/2 
vicinity of the nucleus (see, e.g., [6:]). For instance, 



/pi/2 ('^) 



"1/2 



(R) — Zax 



1 



(18) 



where x = r/ R. 

At r = R we obtain 



/pi/2(0 ~ 9pi,AR) 7;Zax 0.8 



(19) 



and after some transformations we find the correspon- 
dence 



5m.pR ^ 

ZOL 



X 10^ — — Ct 



(20) 



The final coefficient connecting Cp and Ct in Eq. ([^D]) 
differs by a factor of 6/5 from that obtained in ^22,]. 

Given the results obtained in this work in the RPA 
approximation for and using Eq. (I20p we can find the 
values of d^. These values are listed in Table Hill in the 
entry "Rescaling (RPA)" . As is seen from the table there 
is excellent agreement between calculated and rescaled 
values. It is worth noting that Eq. ([20]) turns out to be 
insensitive to Z. For the comparably light Xe and for 
heavy Ra the agreement of the numerical results and the 
results obtained with use of Eq. (PO)) is equally good. It 
means that Eq. (I^Dl) works well for atoms with different 
Z. 



TABLE HI: The values of dft in units (IO'^^Cp (o-iv)|e| cm) 
obtained in DHF and RPA approximations are presented. 
The results are compared with other data. 





i29Xe 


171 Yb 




2"Rn 


Ra 


This work (DHF) 


1.3 


-2.4 


-8.7 


17.3 


-13.0 


This work (RPA) 


1.6 


-11.5 


-21.3 


21.0 


-63.7 


Rescaling (RPA)" 


1.6 


-11.3 


-21.3 


21.3 


-64.7 



"These numbers are obtained using the correspondence 
Cp <-> Ct given by Eq. (|20p as explained in the text. 



TABLE IV: The values of dlC in units 
(lO-^fS/delfm^)] |e|cm) obtained in DHF and RPA 
approximations are presented. The results are compared 
with other data. 





i29Xe 


Yb 


199Hg 


2"Rn 


225 Ra 


This work (DHF) 


0.29 


-0.42 


-1.2 


2.5 


-1.8 


This work (RPA) 


0.38 


-1.9 


-3.0 


3.3 


-8.3 


Ref. [7] (RPA) 


0.38 




-2.8 


3.3 


-8.5 


Ref. [8] (RPA) 




-1.9 








Ref. [23] (RPA) 






-5.07 







2. The Schiff moment 



The contribution to the atomic EDM due to the Schiff 



moment H is naturally to determine as 
df*^(I//), where df*^ is given by 



"at 



jSM 



a.k 



£k 



(21) 



In explicit form the ME of the operator {r p{r)) is given 
in Appendix lAl 

The results for the Schiff moment contribution to the 
atomic EDMs are presented in Table IIVI There is very 
good agreement with previous calculations 0, [!]■ Few 
percent difference for Hg and Ra is within the accuracy 
of the calculations for these atoms. 

A reason of the discrepancy with the result obtained 
in Ref. [2^ is unclear for us. The authors of this work 
state that a possible reason of the difference between 
their result and those obtained in Ref. is due to elec- 
tron correlations which drastically change the final result. 
However, according to our calculations, the correlations 
included in the CI+MBPT approach change the results 
insignificantly, on the level of 10-15% as compared to 
the RPA calculations. The same conclusion was made 
in Ref. ^] . Below we will return to this problem. 



3. The electron EDM 

Usually in experiments searching for atomic EDMs 
atoms are placed in an external electric field Ecxt • It leads 
to appearance of the interaction datEext = — '^e757Eoxt- 
The operator (757) is P-even and T-odd. If the hyperfine 



6 



TABLE V: The values of dft in units (de x 10"") obtained 
in DHF and RPA approximations are presented. The results 
are compared with other data. 





i29Xe 


171 Yb 


199Hg 


2"Rn 


Ra 


This work (DHF) 


0.85 


1.0 


4.9 


-11 


-11 


Ref.[24] (DHF) 


0.86 




5.1 






This work (RPA) 


1.0 


5.1 


12.5 


-13 


-55 


Rcf.f24] (RPA) 


1.05 
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interaction, i?hf, is accounted for, a contribution to the 
atomic EDM caused by this operator appears already in 
the second order of the perturbation theory and looks as 
follows 



2E- 



K 



|gM|g) (^1475710) 
Eq — Ek 



(22) 



But as it was shown in [22], this contribution has the 
magnitude ~ deZ^a^/mp and is negligibly small. We 
disregard it in this work. 

For calculating the contribution of the electron EDM 
(described by the operator Hb) to the atomic EDM we 
use the same approach as for studying the e-N interac- 
tion. We denote the contribution to the atomic EDM 
due to Hb as dfj = d^^{l/I). Using Eq. H]) and replac- 
ing H Hb, after simple transformations we arrive at 
the following expression 



SrririC 



a,k 



(fc|| r||a)(fc||gg||a) 

£k — £a 



(23) 



Here H"^ is the electronic part of the operator Hb- The 
explicit expression for the ME of the operator H'^ is given 
in Appendix [Al 

The results obtained for this contribution to the atomic 
EDMs is presented in TablelV] Comparing our results ob- 
tained for Xe and Hg with those of Martensson-Pendrill 
and Oster we see very good agreement between them. It 
is seen that the RPA corrections change the results ob- 
tained in the DHF approximation for Yb, Hg, and Ra 
very significantly. For Yb and Ra the DHF and RPA re- 
sults differ by a factor of 5. This behavior is quit similar 
to what we found for d^, d^, and dft*^ (see Tables IH 
IHIl and llVp . It is not surprisingly if we take into account 
the similar nature of all these P,T-odd operators. To the 
best of our knowledge there are no other available data 
for Yb, Rn, and Ra. 

There is also the third order contribution to atomic 
EDM proportional de [l^l which is actually larger than 
the second order contribution discussed here. We will 
consider it in a separate publication. 



B. CI+MBPT 

We use the CI-I-MBPT approximation for the calcula- 
tions for Yb, Hg, and Ra atoms, which have two s elec- 



TABLE VI: The removal energies for both 6s electrons for 
Yb and Hg and both 7s electrons for Ra are presented in the 
first row of each respective atom. Energies of excited states 
are presented (in cm"'^) in respect to the ground state. The 
results are obtained in the CI and CI-I-MBPT approximations. 



Yb 



Hg 



Ra 



Config. 


Level 


CI 


CI-fMBPT 


Experime 




'So 


138795 


148707 


148712 


6s6p 


3 r)o 


14368 


17562 


17288 


6s6p 


3po 


15029 


18251 


17992 


6s6p 


3 r)o 


16537 


19995 


19710 


6s6p 


Ipo 


24215 


25715 


25068 


6s2 


'So 


213969 


235409 


235469 


6s6p 




30676 


37537 


37645 


6s6p 


3po 


32446 


39375 


39412 


6s6p 




36541 


44405 


44043 


6s6p 


Ipo 


48195 


54116 


54069 


7s2 


'So 


116301 


124316 


124416 


7s7p 


3 po 

Jo 


10706 


13108 


13078 


7s7p 


3po 


11587 


14003 


13999 


7s7p 




13892 


16693 


16686 


7s7p 


Ipo 


19011 


20597 


20716 



trons above closed shells. For the noble gases such as Xe 
and Rn the RPA is known to be good for describing their 
properties (see, e.g., [1^) and there is no need to use the 
CI-I-MBPT for them. 

We start the discussion of the properties of Yb, Hg, 
and Ra from the results obtained for the low-lying ener- 
gies of these atoms. In Table |VT] we present the energy 
level values obtained in the pure CI and the CI-I-MBPT 
approximation. As is seen from Table IVII the removal 
energies of the two s electrons differ by ~ 10% from the 
experimental values at the CI stage. An accounting for 
the MBPT corrections leads to almost ideal (better than 
0.1%) agreement between these quantities. The energies 
of the excited states (calculated relatively to the ground 
state) are also noticeably improved at the CI-I-MBPT 
stage. The differences between theoretical and exper- 
imental values do not exceed 3%. These results indi- 
cate the accuracy of wave functions produced at different 
stages of the CI-HMBPT method. 

The calculations of the atomic EDM dat in the 
CI-I-MBPT approah is more complicated than in the RPA 
method. Again we start from Eq. ([1]) keeping in mind 
that the summation in this equation is going over many- 
electron states. Following [2^[27l| dat can be divided into 
two parts. 



dat = dl. 



(24) 



where d^t includes excitations of valence electrons and 
dl™° includes excitations of core electrons and a correc- 
tion to da™® that appears because of possible excitations 
of core electrons into the closed valence s shell, what is 
forbidden by the Pauli principle. Note that the correc- 
tion restoring the Pauli principle is not small. In certain 
cases it constitutes ^ 50% of the total core contribution. 
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TABLE VII: The valence and core and contributions to d^t (in 
10-2°CT(fTiv)|e| cm), df^ (in IQ-^^Cp (o-jv>|e| cm), (in 
10""[S/(|e|fm^)] |e|cm), and dg (in 10"*de) are presented. 
The entry "Total' means the sum of the valence and core 
contributions. The results obtained in the RPA method {V^ 
approximation) are given for comparison. We denote this 
entry as "RPA" . 



TABLE VIII: Composition of the nuclear spin. 









dft 




dg 


^''^Yb 

X U 




-4.24 


-14.2 


-2.43 


6.24 




core 


0.54 


1.8 


0.31 


-0.79 




Total 


-3.70 


-12.4 


-2.12 


5.45 




RPA 


-3.37 


-10.9 


-1.95 


5.05 




val. 


-5.71 


-20.5 


-2.95 


12.0 




core 


0.59 


2.1 


0.32 


-1.3 




Total 


-5.12 


-18.4 


-2.63 


10.7 




RPA 


-5.89 


-20.7 


-2.99 


12.3 




val. 


-20.6 


-75.0 


-10.25 


-65.2 




core 


3.0 


10.8 


1.41 


9.5 




Total 


-17.6 


-64.2 


-8.84 


-55.7 




RPA 


-16.7 


-61.0 


-8.27 


-53.3 



With the wave functions obtained from Eq. pS]) . the 
valence part d^^ is c omp uted with the Sternheimer psj 
or Dalgarno-Lewis [29] method implemented in the 
CI-l-MBPT-l-RPA framework. Given the wave function 
|0) and its energy Eq, we find an intermediate-state wave 
function from the inhomogcneous equation 

K 

(r.)efflO). 



1 



HeS — Eq 

Given <i>in we can be compute as 
<t-2(0|(r.)efr|$i, 



(25) 



(26) 



An additional contribution d'^"^ 



coming from particle- 
hole excitations of the core is incorporated in the frame 
of the RPA approach discussed above. 

In Table IVIII we list the results obtained for the va- 
lence and core contributions to the EDM of the atoms. 
For completeness and comparison we also present the val- 
ues obtained in the RPA calculations. As is seen from the 
table the largest differences between the results obtained 
in the RPA and CI+MBPT methods occur for Hg. But 
even in this case these differences do not exceed 15%. 
Taking into account the similar nature of the P,T-odd 
operators considered in this work we can expect that the 
relative difference between the results found in both ap- 
proximations would be approximately the same for all 
operators for a given atom. As we see from Table IVIII 
this condition is fulfilled. This is a good consistency test 
of the calculations. 

In this way we arrive at the conclusion that the elec- 
tron correlations do not affect the final results too much. 
Based on this observation we estimate the accuracy of the 



Nucleus 


Neutron state 




{^'"^} 




^^«Xe 
«9Hg 


Pl/2 


1 

-1/3 


0.76 
-0.31 


0.24 
-0.03 



values "Total" listed in Table IVH is at the level of 15- 
20%. Similar accuracy is expected for Xe and Rn atoms. 



V. SHELL MODEL OF THE NUCLEUS 

In certain cases it is possible to impose additionally the 
constraints on the strength of the couplings for protons 
and neutrons. Below we will obtain such constraints for a 
simple case of the spherical shell model of a nucleus. This 
model describes well the nuclei of ^^^Xe and ^^^Hg. The 
former contains a neutron in the Si/2 state and the lat- 
ter contains a neutron in the pi /2 state above the closed 
shells. Although the nuclei in the atoms have a valence 
neutron it is possible to deduce also the proton contribu- 
tion to the total nuclear spin using the information on the 
nuclear magnetic moments. If we assume that the mag- 
netic moment of the nucleus is composed entirely from 
the spin magnetic moment of the valence neutron and 
spin magnetism of polarized nuclear core, then 



M = Mn('Ti")>+Mp(a(?')) 



(27) 



Here fip « 2.793 and w —1.913 are the magentic 
moments of the proton and neutron expressed in nuclear 
magnetons. Numerical estimates show that the main con- 
tribution to the nuclear magnetic moments /i of ^^^Xe 
and ^^^Hg comes from the neutron and proton spin con- 
tributions. This is due to that the neutron orbital con- 
tribution is zero and the proton orbital contribution (for 
the orbitals with a low orbital momentum) is small in 
comparison with its spin contribution. Neglecting the 
spin-orbit interaction leads to conservation of the total 
spin which is equal to the average spin of the neutron 
above the unpolarized core {a^^^). Taking into account 
that in the spherical shell model the nuclear spin I is de- 
termined by the total momentum of the unpaired nucleon 
we can write 



with 



.(o)\ 



1, 

-Ilil- 



(28) 



/ = ;/ + l/2 
1), / = Z/-l/2 



where li is the orbital quantum number of the valence 
nucleon. Using Eq. P5|) . we determine (cri"'') and (cri^^) 
for observationally relevant cases of ^^^Xe and ^^^Hg as 
shown in Table EID 
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As is seen from Table IVllIl the contribution of the pro- 
ton spin into the total nuclear spin of ^^^Xe is as high 
as 30%, and therefore the proton couplings Cj. and Cp 
(see the equations below Eq. ^) are also limited in the 
experiments searching for the EDM. For ^^^Hg the limit 
on the proton couplings is 10 times weaker than the limit 
on the neutron couplings. 

Note that the simple shell model of a nucleus consid- 
ered above is hardly applicable to the nuclei of ^^^Yb, 
^^^Rn, and ^^^Ra. The problem is that the nucleus of 
^^^Yb is quadrupole-deformed and the nuclei of ^^^Rn 
and ^^^Ra are octupole-deformed. For this reason more 
sophisticated nuclear models are required for a proper 
description of these nuclei. 



VI. CONCLUSION 

In conclusion, we have carried out calculations of dif- 
ferent contributions to the EDMs of ^^SXe, "^Yb, ^^SHg, 
2"Rn, and 225Ra by the RPA and the CI-hMBPT meth- 
ods. Two of these contributions are due to tensor- 
pseudotensor and scalar-pseudoscalar e-iV P,T-odd inter- 
actions, and two more contributions are caused by the 
nuclear Schiff moment and the electron EDM. The rec- 
ommended values for noble gases Xn and Ra are based 
on the results obtained in the RPA calculations while for 
Yb, Hg, and Ra they based on the calculations carried 
out in the frame of CI+MBPT-fRPA approach. These 
numbers are gathered in Table ITXl 

Finally, using the results obtained in this work for 
atomic mercury and recently obtained new upper bound 
Mat(^^^g)| < 3.1 X 10-29 |e|cm Q we are able to find 
limits on CP violating parameters Ct, Cp, S. We do not 
put a constraint on the coupling constant de because the 
contribution of the electric dipole electron moment to the 
atomic EDM in the third order of the perturbation theory 
is known to be ~ 10 times larger than its contribution in 
the second order . It means that a dominating contri- 
bution from the electron EDM still has to be considered. 
This is a subject of another work. Applying again the 
spherical shell model to ^^^Hg and, respectively, having 
((Tjv) = —1/3(1//) we arrive at the numbers listed in 
Table H 

We would like to thank M. Kozlov for useful remarks. 
This work was supported by Australian Research Coun- 
cil. S.G.P. was also supported in part by the Russian 
Foundation for Basic Research under Grants No. 07-02- 
00210-a and No. 08-02-00460-a. 



APPENDIX A 

To calculate the MEs of the P,T-operators we define 
the one-electron wave function \a) = ■(/'a(r) as follows 



V'Q(r) 



fair) n^j^rnA^) 

iga{r) (n) 



(Al) 



where la = 2ja - la- 

The MEs of the P,T-odd operators characterizing the 
tensor-pseudotensor and pseudoscalar-scalar interaction 
are given by the following expressions 



nbKbWil p{r)\\naKa) = (^bl |Ci | l^a) X (A2) 

{fbgaiKa- Kb+ 1) +fagb{Kb- K.a+ I)} pir^"^ dr, 



{nbKb\\n {dp/dr)jQ\\naK.a) = {Kb\\Ci\\Ka) 



X J {9b9a - fbfa) ^dr. 



(A3) 



If the nuclear density distribution is given by Eq. ^ (see 
the main text) Eq. (jA3p is further simplified leading to 



{nbKb\\n{dp/dr)-fo\\naKa) = 

/ 11/^11 \ ^ ( 9b9a fbfa 



(A4) 



r=R 



where spherical harmonics Cim arc defined as 



C;m(n) = 



47r 
2/ + 1 



(A5) 



K — {I — j)(2j + 1) and the reduced ME (K6||Ci||Ka) is 
given by 

(nbWCilK) ^ ah + la + I) y< (A6) 

3b ja 1 ' 



(-l)^''+^/V(2ja + l)(2j5 + l)(^_f/2 {J^ J 
1, if a; is even 



where f(x, — >, n -r • jj 
^ ' 0, if X IS odd 

The reduced ME characterizing the nuclear Schiff mo- 
ment is given by 

{ubUbWr p{r)\\naKa) = (Kb||C'i||KQ) X (A7) 
{9b9a +fbfa}r^ p{r)dr. 



The reduced ME of the operator Hg can be repre- 
sented as 



{nbKb\\H'l^\\naKa) = ~{Kb\\Ci\\Ka) X 



(A8) 



' dr 

{fbga{Kb- K.a+ 2) + faQbi^a^ Kb+ 2)} 

r 

/■^ r^dr 

2 J {fb9a{nb~ Ka~ ^) + fa9b{l^a- Kfc - 1)} 
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TABLE IX: The recommended values of the contributions to atomic EDM. 





i29Xe 


Yb 




2"Rn 


Ra 


d^t(10-""CT(o-Ar)|e| cm) 
dg{10-''^ Cp {a N)\e\ cm) 
dl^{10-^^[S/{\e\im^)]\e\cm) 
dft(10'*de) 


0.57 
1.6 

0.38 
1.0 


-3.7 
-12 
-2.1 
5.5 


-5.1 
-18 
-2.6 

11 


5.6 
21 
3.3 

-13 


-18 
-64 
-8.8 

-56 



TABLE X: Limits on CP violating parameters Ct, Cp, and 
S based on new experimental limit for |dat(^^^g)l < 3.1 x 
10-29 |e| cm 3 



Parameter 


Limit 


Ct 


1.9 xlQ-" 


Cp 


5.2 xlO"'^ 


5(le|fm3) 


1.2 X 10-12 
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